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Abstract 

In this paper, we propose a realistic mathematical model taking into account 
the mutual interference among the interacting populations. This model at¬ 
tempts to describe the control (vaccination) function as a function of the 
number of infective which is an improvement over the existing susceptible- 
infective epidemic models. Regarding the growth of the epidemic as a nonlin¬ 
ear phenomenon we have developed a neural network architecture to estimate 
the vital parameters associated with this model. This architecture is based 
on a recently developed new class of neural networks known as co-operative 
and supportive neural networks and it involves a preprocessing of the input 
data and this renders an efficient estimation of the rate of spread of the epi¬ 
demic. It is observed that the proposed new neural network outperforms a 
simple feed forward neural network and polynomial regression. 
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1. Introduction 

Simple epidemic models describe the spread of an infectious disease among 
individuals within the population. After a period, two groups are formed in 
the population: those who have not acquired the disease but are likely to 
contract (susceptible population) and those who are infected (infective pop¬ 
ulation), capable of spreading the disease. The fundamental characteristic of 
the epidemic is that susceptible individuals contract the disease only by get¬ 
ting in contact with the infective individuals. Also the average time constant 
or the latency time of the disease depends on the nature of the epidemic. 
Further, cured individuals do not contract the disease again during the same 
period. Generally, such epidemics will be treated by appropriate vaccina¬ 
tion and/or other efforts, which may be viewed as control efforts to contain 
the spread of disease. The vaccination effort is regarded as a parameter in 
the mathematical models. The estimation of the rate at which susceptible 
individuals become invectives is generally a difficult question for the held 
scientists engaged in this activity. So it is desirable to have a realistic math¬ 
ematical model that describes the dynamical interactions between these two 
classes of populations, such as the susceptible and the infective. For a few 
earlier studies on epidemiological problems we refer the readers to (121 Eli, 
[3 OSl nsi ESI DS]). This is the starting point for our investigations in this 
paper and accordingly, we propose a mathematical model to describe the 
simple dynamics of the interacting groups of the populations. 

We introduce the nonlinearities in the interacting populations through 
mutual interference parameters. Our main result provides suitable ranges 
for these parameters. Regarding the growth of the epidemic as a nonlinear 
phenomenon, we have developed a neural net work architecture to estimate 
the vital parameters of the model. Convergence of the neural network train¬ 
ing is an important problem in case of excess data samples being available. 
We propose a new methodology to train a neural network for data intensive 
applications by employing clustering techniques. 

The present paper is organized as follows. In Section 2, we describe our 
mathematical model. In Section 3, we derive conditions that ensure the 
existence and uniqueness of continuable solutions for the model equations. 
A question of importance to the held scientists is the determination of the 
rate of spread of the epidemic and we utilize a recently developed neural 
network architecture [17] to estimate this rate. A related algorithm and the 
new neural network architecture with K-means clustering are discussed in 
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Section 4. Simulation results are presented in Section 5. The rate of spread 
of the epidemic is determined using the neural network architecture developed 
in Section 4 and these results form the content of Section 6 . Conclusions and 
discussion follow in Section 7. 

2. The Model 

Our model is based on the consideration that the process of the spread 
of the susceptible population turning into infective population is a nonlinear 
phenomenon. The following are the underlying biological principles. 

1. The total population is hxed and initially every individual is susceptible 
to the disease. 

2. The disease is spread through the direct contact of susceptible individ¬ 
uals with the infective individuals. 

3. Every individual who has contracted the disease and has recovered is 
regarded as immune. 

These principles when translated into the mathematical framework yield 
the following system of ordinary differential equations and these equations 
describe the dynamics of the interacting populations, 

- S{xi) ( 2 . 1 ) 

- 7^(a:2) 

where , =^, Xi represents the number of susceptible individuals, X 2 rep¬ 
resents the number of infective individuals, (5 is the infection parameter, 7 
denotes a parameter related to the average time constant of the disease. The 
function S describes the control input and it is assumed to be proportional to 
the vaccination effort. Also, the function P corresponds to those individuals 
who have contracted the disease and recovered (regarded as individuals with 
acquired immunity). Further the functions S and P are assumed to satisfy 
the following mathematical conditions: 

S(0)>0.A>0;P(0) = 0.^>0 (2.2) 

The parameters mi and m 2 appearing in (2.1) represent the indexes of the 
interacting populations arising out of the non-linear considerations of the 
epidemic phenomenon. Prototypes of the model (2.1), in which the nonlinear 
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interactions mi = m 2 = 0 and when a fixed (constant) control effort is 
invoked, reduce the above model equations to 

Xi = —(3xiX 2 — u (2.3) 

X2 = I3xiX 2 - 'yx2 

and this model has been studied in ([I]). More recent studies on the better 
vaccination efforts may be found in the recent papers ([121 na dsi [n]). 

Clearly the present model (2.1) is more realistic as the control effort es¬ 
sentially varies with regard to the size of susceptible population rather than 
being fixed. The model (2.3) rests on the simple considerations that if each 
infected individual converts one susceptible into infective, then the total num¬ 
ber of susceptible population converted into infective population would be 
and this describes the interactions between these populations (simple nonlin¬ 
ear interactions). Our model instead of considering simple nonlinear interac¬ 
tions also addresses the sub linear interactions between the two populations. 

3. Existence and Uniqneness of Continuable Solutions 

In this section, we propose to consider the model equations (2.1) and 
examine the qualitative properties of its solutions. From the biological point 
of view, the system (2.1) describes the dynamical interactions among the two 
classes of populations such as the susceptible and the infective populations. 
Clearly, the qualitative study of solutions of this system depends on ensuring 
conditions that are sufficient to guarantee the existence of solutions for initial 
value problems associated with (2.1). Usually this existence of solutions 
is determined in a finite interval and the solutions are continued on their 
maximal intervals. This approach yields continuable solutions. Often the 
inherent dynamics of the system requires one to pick up a specific point on 
the trajectory and to move in both forward and backward directions. From 
this discussion, a mathematical treatment of the model equations requires one 
to obtain conditions for the existence and uniqueness of continuable solutions 
for the system (2.1). An analogy with ecological problems would render one 
to regard the community of infective individuals as a sub-population affecting 
the survival of the susceptible individuals. One, aspect of the dynamics of 
community interactions would be mutual interference among the interacting 
sub-populations, which in our model (2.1) is represented by the parameters 
mi and m 2 . It is known that mutual interference is a stabilizing process 
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mu fT0]t A question of interest is to describe conditions leading to the 
persistence/survivability of interacting populations. We observe that the 
mutual interference introduces sub-linearities into the system (2.1) and it 
is known that initial value problems for systems with sub-linearities have 
continuable solutions but these solutions are not unique. (naci), and hence 
may not be regarded as a dynamical system (BE!). Our hrst result in 
this direction is to hnd suitable ranges for parameters rrii and m 2 so that the 
solutions of (2.1) form a dynamical system in the sense described above 

Theorem 3.1. Consider the system of equations given by 


Xi = gi{xi,X2), (3.1) 

where (rj(0) > 0, The functions gi : R~^ —)■ = [0, cx)) are continuous, 

fori = 1,2, that is g^ G C{R~^). 

Proof. Assume that the following conditions are satished: 

(HI) There exists constants > 0 such that hj G C{R^) where hj{xi, X 2 ) = 

(H2) iP jf; 


Then the solutions of the system (3.1) form a dynamical system in the 
sense of (i)- It is easy to see that the proof is a slight modihcation of a 
result of (|1]). We now apply the content of Theorem 3.1 to our model (2.1). 
Consider the following transformation of the variables for the system (2.1) 
given by 


Ml 

U2 


X 


1—mi 
1 


X 


l-m2 

2 


(3.2) 


Xi 


X2 


1 


u 


1 —m,]^ 

1 


l-’"2 

^2 





1 — m 

Xi 

- Ui 


m 2 

7719 

1 — m 

X 2 

- U 2 


1 


2 


(3.3) 
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Differentiating equation (3.2) 

til = (1 — rni)x^^xi Xi = 

ti 2 = (1 - m2)x^^X2 X2 = 


1 — mi 
1 

1 - m2 


x^tii 


X'^'^U2 


(3.4) 


Substituting (3.4) in (2.1) and after some simplifications, we get 

r ^2 -mi 1 ■ 

■ui = —(1 —mi) /3 M2 ^“"*2 + tip-™! S'(tip-’"i) 

r rni —m2 1 

■ U 2 = (1 — m2 ) /3 Ui — 7 M 2 ^“’"2 P { U 2 ^“™2 ) 


(3.5) 


In order to verify the hypotheses (HI) of Theorem 3.1, we need to show 


that 

exists, in which 

hi{Xi,X2) = Xi~'^^ 


lim hj{xi,X 2 ) 


4>0+ 


(3 xi^^X2^‘^ + *S'(a;i) = (3 0 : 2 ™^ + xi S{xi) 


h 2 {Xi,X 2 ) = X 2 


13 Xi^^ X2^^ — 7 P{X2) 


[3 Xi^^ — 7 3(2 P{X2) 


From the assumptions (2.2) on the functions S and P it is clear that the 
hypothesis (HI) of Theorem 3.1 is verihed. Also, it follows easily that the 
hypotheses (H2) of Theorem 3.1 are verihed provided the parameters m^ 
satishes the inequalities 

2mfc — 1 > 0 


□ 


Finally, an application of Theorem 3.1 yields the conclusion that solutions 
of the model (2.1) describe a dynamical system provided 

1 

mfc > - 

. Henceforth, we designate the inequalities 

1 

mfc > - 

as admissible values of the parameters mi and m 2 . 
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4. Neural Network Architecture 

As mentioned above, a question of practical interest has been to determine 
the rate at which the susceptible become infective. In this section a new 
neural network architecture with back-propagation algorithm is proposed. 

The training data is obtained by solving the model equations (2.1) using 
Matlab. A simple neural network without a proper preprocessing will not 
converge due to the non homogeneity in the data sets. Therefore a prepro¬ 
cessing step with K-means clustering is introduced to create homogeneous 
data sets before training the network. The neural network is trained using 
Xlminer Microsoft Excel 2003 plug-in. 

4 . 1 . Methodology 

The k-means clustering works on the expectation of maximization algo¬ 
rithm to hnd the centers of natural clusters in the data. It assumes that 
the object attributes form a vector space. The objective is to minimize total 
intra-cluster variance, or, the squared error function 

k 

V = 

*=1 XjGSi 

where there are k clusters Si,i = 1, 2 ,..., fc and fii is the centroid or mean 
point of all the points Xj G S^. The initial number of clusters is specihed as 
three based on the visual inspection of the data sets, for generating clusters 
from the data using k-means clustering. The following Figure 1 illustrates 
the operating mechanism for training a neural network. The neural network 
architecture is depicted in Figure 2. 


Figure 1: Operating Mechanism forof the proposed Learning paradigm for cooperative 
networks. 
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Figure 2: Cooperative and supportive multi layer feed forward neural network architecture 
for computing susceptible for a given infective population. 


5>nscei>lihles 



Tnfec lives 


K-means clustering is a parametric technique where it is required to pro¬ 
vide the number of clusters as a parameter. To know how well-separated the 
resulting clusters are, a Silhouette plot is constructed using cluster indices 
output from k-means. 

The Silhouette plot displays a measure of how close each point in one clus¬ 
ter is to points in the neighbouring clusters. This measure ranges from -(-1, 
indicating points that are very distant from neighbouring clusters, through 
0 , indicating points that are not distinctly in one cluster or another, to -1, 
indicating points that are probably assigned to the wrong cluster. From the 
Figure 3 it is derived that the optimal number of clusters is three as silhouette 
plots of four and hve clusters shows very low Silhouette values. 

4-2. Initial Training Procedure 

Initially the entire susceptible and infective populations are standardized 
to mean zero and standard deviation one. Using k-means clustering algorithm 
the data sets are clustered. The set of susceptible and infective popnlations 
in each cluster is then partitioned in to training data and testing data. The 
training data is given as input to the neural network. When the training is 
complete the testing data is given as inpnt and outputs are computed. 
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Figure 3: Silhouette plot with five, four and three clusters generated using k-means clus¬ 
tering. 



The three layer neural network architecture with one input, hve hidden 
neurons and one output is considered for training. Each layer has one bias 
node. The tan hyperbolic function is used as a threshold function for each 
of the neurons. The network is trained using the back propagation learning 
algorithm with momentum, which has been found effective. The procedure 
and the network architecture designed are shown in Figure 4. The data 
used for training is generated from the mathematical model (2.1) through 
numerical solutions using Matlab. 

The network is trained to predict the number of susceptible population 
given the number of infective population to the input neuron. The purpose of 
predicting susceptible is to estimate the rate of spread of the epidemic. Dur¬ 
ing the training phase the error generated due to differences in the predicted 
and actual susceptible populations is propagated in the backward direction 
and weights are adjusted. The network is trained once the mean square error 
reaches a dehned required value. When the training is complete the infective 
in the testing data set is given as input and the susceptible are computed. 
Figure 5 depicts the structure of the trained net work which computes the 
susceptible population with infective population as input. 

The training and testing are done for each of the clusters and the suscep¬ 
tible are computed for the test data sets for all the clusters. 
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Figure 4: Neural network architecture for training. 

SLjeeprb« ^arr 

nng data 



Figure 5: Testing procedure for the estimation of the susceptible given the infective as in 
input. 


from test 
dataset 


1 



■ 0 " 


l-tOT)' 


5. Simulations 

The first part of this section deals with the results of the model (2.1) 
simulated for various admissible values of the parameters mi and m 2 . Model 
1 considers a situation in which sub linearity between the susceptible and 
infective populations is not taken into account and with hxed vaccination 
effort. 

5.1. Model 1 

This example deals with the case in which mi = m 2 = 1, m = 10, f3 = 
0.0001, 7 = 0.8, Susceptible population = 10,000, and infective population 

= 10. This corresponds to the model studied in (pQ) 

Xi = —f 3 X1X2 — u, 

X 2 = /3 Xi X 2 - ^X2 


Figure 6 shows the interactions between the susceptible and the infective 
populations generated using Matlab software. 

5.2. Comparison of Neural Network and Regression Analysis for Model 1 
We have applied statistical methods such as regression analysis for Model 
1 and compared the performance of the neural network with the statistical 
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Figure 6: Interactions between susceptible and infective population using Model 1. 


SUSCEPTIBLES Vs INFECTIVES OF MODEL1 SOLUTIONS 



methods. Neural Network architecture with conhguration 1-5-1 is designed. 
When k-means clustering is applied to the input data three clusters are ob¬ 
tained. Each individually clustered data is given as input to the three neural 
networks for training using back-propagation with momentum learning al¬ 
gorithm. Also a linear and quadratic regression analysis is carried and the 
results are given in Figure 7. 


Figure 7: Susceptible individuals computed from neural network and regression analysis 
giving infective as inputs for Model 1. 



Figure 7 shows the actual susceptible and the susceptible estimated from 
the neural network and regression analysis. It may be seen that the estima¬ 
tion by the neural network is better when compared with that obtained by 
using regression method. 
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5.3. Model 2 


xi = —jd X2^'^ — S{xi) 

X2 = (d Xi^^ X2"^^ — 3 P{x 2 ) 

This model is an improvement over the Model 1 in the sense that the vaccina¬ 
tion effort is based on the infection rate and is dynamic in character. Also, 
the interactions between the populations are not necessarily linear. This 
model has the following parameters 


mi = 0 . 8 , 

m 2 = 0.7, 

s{x) = {xi^’'^ / {vaccination ef fort + Xi^'"^)) 

P{x) = 

fd = 0.01, 7 = 0.04, initial susceptible = 1000, and initial infective = 10. 

The interactions between the susceptible and infective populations are 
shown in the Figure 8 which has been generated using Matlab software. The 
vaccination effort has rendered decline in the infective population well ahead, 
as observed in Model 1. This clearly establishes that the vaccination effort 
given by Model 2 is more realistic. 


Figure 8: Interactions between susceptible and infective individuals during the spread of 
the epidemic using numerical simulation in Matlab for Model 2. 



TIME 
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5.4- Neural Network Training with Infective as Input and Susceptible as ex¬ 
pected Output for Model 2, 

Neural network architecture similar to the architecture discussed in Sec¬ 
tion 5.2 is considered for training the network. Linear and quadratic regres¬ 
sion analyses are carried out on this data and the results are plotted in the 
Figure 9. From the Figure 9 it may be observed that the neural network esti¬ 
mation of susceptible is better than those obtained by the regression method. 


Figure 9: Susceptible individuals computed from neural network and regression analysis 
with infective population as input for Model 2. 



TIME 


6. Estimation of the rate of Spread 

In the section, we employ the neural network architecture and estimate 
the rate of spread of the epidemic and compare the same with the actual 
calculated rate. The rate is the product 

/3 


6.1. Model 1 

This example shows the rate calculated from a neural network and the 
actual calculated rate ( 
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) are in good agreement than that obtained from polynomial regression. Sim¬ 
ulations are conducted with the following parameter values mi = m 2 = 
1 , /? = 0 . 001 . 

The Figure 10 shows that the rate estimated by the neural network is 
very close to the actual rate compared to the one estimated by using the 
regression methods. 


Figure 10: Estimated rate from neural network, regression analysis and actual calculated 
rate of spread of epidemics Model 1. 



6.2. Model 2 

In this example Model 2 is considered with mutual interference param¬ 
eters and the rate of spread of the epidemic is calculated from the neural 
network and also the polynomial regression method for the following values 
of the parameters mi = 0.8, m 2 = 0.7, fd = 0.01. 

In all these examples it is observed that the results obtained through the 
neural network are closer to the actual results than those obtained by using 
the polynomial regression analysis. 

7. Conclusions and Discussion 

This work attempts to estimate the rate of spread of an epidemic un¬ 
der realistic conditions, which is a problem of immense concern for the held 
scientists. It has been realized that a realistic mathematical model is es¬ 
sential to deal with this problem. Accordingly, a mathematical model has 
been proposed and certain mathematical questions such as the existence and 
uniqueness of continuable solutions to the model equations (which ensure 
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Figure 11: Estimated rate from neural network, regression analysis and actual calculated 
rate of spread of epidemics. 



that the system describes a dynamical system) have been discnssed. Follow¬ 
ing a new learning paradigm using k-means and cooperative neural networks 
the rate of the spread of epidemic has been estimated. Results of earlier 
work have been derived from the present work. Simulation results exhibit 
the decline in both susceptible and infective populations with increased con¬ 
trol effort; there by implying that the populations contracting the epidemic 
are cured and/or the vaccination effort makes them immune. A non-linear 
regression analysis is also carried out for all the models and it is concluded 
that the non-linearity in the data is better adapted by the neural network 
approach than the regression analysis. The neural network out performs the 
regression analysis by predicting the rate very close to the actual rate. It 
is hoped that this work paves way for better understanding of the simple 
epidemic phenomenon. 
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